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Abstract 

In this study we combine the multiwavelength ultraviolet - optical (Solar Dy¬ 
namics Observatory, SDO) and radio (Nobeyama Radiohcliograph, NoRH) 
observations to get further insight into space-frequency distribution of oscil¬ 
lations at different atmospheric levels of the Sun. We processed the observa¬ 
tional data on NOAA 11711 active region and found oscillations propagating 
from the photospheric level through the transition region upward into the 
corona. The power maps of low-frequency (1-2 mHz) oscillations reproduce 
well the fan-like coronal structures visible in the Fe IX 171A line. High fre¬ 
quency oscillations (5-7 mHz) propagate along the vertical magnetic held 
lines and concentrate inside small-scale elements in the umbra and at the 
umbra-penumbra boundary. We investigated the dependence of the domi¬ 
nant oscillation frequency upon the distance from the sunspot barycentre to 
estimate inclination of magnetic tubes in higher levels of sunspots where it 
cannot be measured directly, and found that this angle is close to 40° above 
the umbra boundaries in the transition region. 
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1. Introduction 


The waves in solar magnetic structures play an important role in pro¬ 
cesses of the energy exchange. Moreover, the properties of the medium, i.e. 
the solar atmosphere, can be determined by studying the wave propagation. 
This affects such problems as the origin of the solar magnetic fields, coro¬ 
nal heating, flares, wave propagation in plasma with different magnetic field 
topology, and others. From the observation driven point of view the prob¬ 
lems above depend on our knowledge of the magnetic held distribution with 
height, plasma velocity, temperature and their variations. Because direct 
measurement of these parameters is a considerable challenge, in practice one 
has to use different tools to get a comprehensive picture of the atmosphere 
structure in active regions and processes therein. The seismological study 
of the active regions is one of such unique tools (Bogdan and Judge, 2006; 
De Moortel and Nakariakov, 2012). Different methods were developed to 
explore magnetic geometry of the sunspots and to estimate the magnetic 
held inclination angle in the transition region and in the corona (Jess et ah, 
2013; Yuan et ah, 2014b; Kobanov et ah, 2015). Coronal seismology pro¬ 
vides a capability to determine the magnetic held strength in coronal closed 
magnetic structures (Nakariakov and Ofman, 2001). Another opportunity is 
radio magnetography (Gelfreikh, 2004). The radio emission is naturally con¬ 
nected with magnetic held providing an instrument to study the 3D-structure 
of the sunspot atmosphere. Joint observations in radio and optical spectral 
ranges have a great potential (Kislyakov et ah, 2006), that has not been 
implemented yet. One can trace the wave train propagation and utilise mag¬ 
netic held extrapolation to the corona to estimate the effective height of the 
radio emission, which can be used to qualify the formation heights of the 
lines in UV and optical spectral ranges. New facilities like the Siberian Solar 
Radio Telescope (SSRT, Lesovoi et ah 2012), the Chinese Spectral Radiohe¬ 
liograph (CSRH, Yan et ah 2013), and others will provide new possibilities 
for observing the Sun’s radio emission in the multifrequency mode. 

Sunspot oscillations were extensively studied in the past decades (see re¬ 
views by Lites (1992); Solanki (2003); Bogdan and Judge (2006); Jess et ah 
(2015)). Recent studies are based on multiwavelength observations and focus 
on determining the wave propagation speed (Felipe et ah, 2010; Abramov- 
Maximov et ah, 2011; Kiddie et ah, 2012; Chandra et ah, 2015), wave dy¬ 
namics (Sycli et ah, 2012; Sycli and Nakariakov, 2014; Zhugzhda and Sych, 
2014), and related phenomena like jets, penumbral filaments, umbral flashes 
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Figure 1: Sunspot NOAA11711. Left panel: SDO/HMI continuum image. Middle panel: 
the sunspot in the SDO/AIA 1700 A band pass. Right panel: magnetic field normal projec¬ 
tion derived from the SDO/HMI vector data. The inner and outer penumbra boundaries 
are outlined by black and white solid lines. The dashed isoline marks the regions where 
the magnetic field inclines at 45 to the surface normal. 


(Yurchyshyn et al., 2014, 2015; Esteban Pozuelo et al., 2015; Madsen et al., 
2015; Louis et al., 2014; Maurya et al., 2013), light bridges (Yuan et al., 
2014a; Yang et al., 2015) and also flares (Sych et al., 2014). Due to the 
circular symmetry, sunspots provide an opportunity to study the wave prop¬ 
agation processes in the presence of magnetic fields. In sunspot umbra, five- 
minute oscillations dominate at the photospheric level, and three-minute ones 
prevail in the chromosphere (Kobanov, 1990; Kentischer and Mattig, 1995; 
Sigwarth and Mattig, 1997; Rouppe van der Voort et al., 2003; Kobanov 
et ah, 2009; Reznikova and Shibasaki, 2012). Oscillations with longer periods 
(10-15 minutes) manifest themselves at the peripheral part of the sunspots. 
Other magnetic field regions, e.g. facula show very different behaviour of the 
spatial-frequency dependence (Kobanov and Pulyaev, 2007; Centeno et ah, 
2009). Detecting and studing the coronal counterparts of these waves and 
quasi-periodic pulsations are relevant for the problem of energy flux transfer 
to the upper atmosphere. 

We continue studing the space-frequency stratification of the oscillation 
processes in the sunspot atmosphere (Kobanov et ah, 2015). We extend the 
analysis with new sunspot data, including those obtained in radio emission, 
and compare the results with magnetic field extrapolation. 
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2. Observational Data and Methods 

The paper is result of the investigation into the sunspot oscillations ob¬ 
served at different wavelengths representing several heights of the solar at¬ 
mosphere from the photosphere to the corona. We used the data from two 
observatories in this study: Solar Dynamics Observatory (SDO, Lemen et al. 
(2012); Scherrer et al. (2012); Woods et al. (2012)) and Nobeyama Radio- 
heliograph (NoRH, Nakajima et al. 1994). The data were selected to cover 
height scale from the photosphere to the corona. 

The sunspot had regular circular shape (Figure 1, left panel). The obser¬ 
vation was performed on 2013 April 6 00:00-04:20 UT and there was no flare 
activity at this period; a Cl.7 flare occurred on 10:58 UT. 

We chose four spectral bands provided by the Atmospheric Imaging As¬ 
sembly (AIA) for the analysis: Feix 171 A and Fexn, xxiv 193 A (coronal 
heights), Hell 304 A (transition region) and 1700 A (upper photosphere, Fig¬ 
ure 1, middle panel). The Hclioseismic and Magnetic Imager (HMI) data 
were used to obtain the magnetic field at the photospheric level in the Fel 
6173 A line, which is formed at the 100-150 km height (Fleck et ah, 2011 ; 
Parnell and Beckers, 1969). The vector magnetic-field components were com¬ 
puted using the code of Very Fast Inversion of the Stokes Vector (Borrero 
et ah, 2011). The SDO data sets were prepared using the aia_prep and 
hrai_prep SolarSoft routines. The sunspot was observed near the central 
meridian (W03S17), which allowed us to minimize the projection effect on 
the analysis of oscillations at different heights. We corrected the magnetic 
field angle to the solar normal within the planar approximation using equa¬ 
tion (Gary and Hagyard, 1990): 

(p = cos'ycostp' + sin'y sirup' cosO ( 1 ) 

where 7 is the angle between the line of sight and the solar normal; tp 1 is 
the observed field inclination angle; 6 is the azimuthal angle corrected for 
the line-of-sight deviation from the central-meridian plane (Figure 1, right 
panel). 

The radio data are the sequence of the 17 GHz NoRH images with the 
20-second effective cadence and the 5" spatial resolution. The images were 
synthesised with the “Koshix” algorithm on the Solar Data Analysis System. 
To increase the signal-to-noise ratio, the number of integration frames for syn¬ 
thesis was set to 10. The images were cropped to include the sunspot area 
(400 x 400"). Then, the images were co-aligned using the cross-correlation 
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technique and smoothed over spatial dimensions. The maximum brightness 
temperature in the sunspot does not exceed 29000 K and the degree of polar¬ 
isation is about 25%. Thus, the polarised emission at 17 GHz in this image is 
most likely due to optically thin thermal bremsstrahlung. The line-of-sight 
magnetic held was estimated by the method proposed by Gelfreikh (2004). 

To study space-height properties of the sunspot oscillations we used the 
method described in Kobanov et al. (2015) that is suitable for circular sunspots 
To calculate the cut-off frequency, we used the equation (Bel and Leroy, 1977; 
McIntosh and Jefferies, 2006; Botha et ah, 2011): 

= 

4:TTV s 

where f c is the cut-off frequency; go=274 ms -2 is the gravitational constant; 
v s is the speed of sound; 7 = | is the adiabatic index; <p(r) is the dependence 
of inclination angle on the barycenter distance. Using Equation 2 and em¬ 
pirical ratio f c =0.81F(r) mHz— 1 mHz for the 304 A line (Yuan et ah, 2014b) 
we estimate the magnetic-field inclination angle [ 79 ] (Kobanov et al. (2015)): 


<p = arccos 


/ 4irv s F(r ) \ 
V 9o1 ) 


( 3 ) 


The oscillations’ power spectrum was calculated using the Fast Fourier 
Transform (FFT). To remove trends from the data we subtracted the time 
series smoothed over a 25-minute time interval from the original series. The 
data set temporal resolution varies from 12 to 45 s for different instruments 
and bandpasses, which allows us to analyse frequencies up to 11 mHz - 
the Nyquist frequency for the HMI data; the Nyquist frequencies for the 
AIA data are 21 mHz (the 1700 A band) and 42 mHz (the 304 A and 171 A 
bands), and that for the NoRH data is 25mHz. To analyse the spatial 
distribution of the dominant frequencies, we obtain FFT power maps in the 
following way. FFT power spectra were calculated for every pixel (0.6 x 0.6 
arcsec 2 ) of the SDO/AIA/HMI intensity images. The calculated Fourier 
spectra were smoothed by convolution with a rectangular window 1 mHz 
wide. The frequency corresponding to the power spectrum maximum value 
was taken as the dominant frequency for the particular spatial domain. To 
analyse wave trains and to measure their propagation speeds the frequency 
filtration was performed by the sixth-order Morlet wavelet. We used the 
algorithm described by Torrence and Compo (1998) to derive the FFT power 
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Figure 2: Spatial distributions of dominant frequencies from the photosphere to the corona. 
The black solid lines mark the inner and outer penumbral boundaries. 


spectra, to wavelet-filter the data and estimate the statistical significance. 
The normalisation by 1 /cr 2 gives a measure of the power relative to the white 
noise level: 3cr 2 correspond to the 95% confidence level, and 4.6cr 2 — 99% 
confidence level (a 2 is the variance of a time series). 

The coronal magnetic held was calculated in the potential and nonlinear 
force-free (NLFF) approach, using SDO/HMI vector magnetograms as input 
data. To remove 180° ambiguity in the direction of hied component transver¬ 
sal to the line-of-sight we used the SFQ method by Rudenko and Anhnogen- 
tov (2014). The potential held was calculated using the method, described in 
Alissandrakis (1981) by applying FFT. The NLFF reconstruction was per¬ 
formed by using the optimization method, proposed by Wheatland et al. 
(2000). This method implementation by Rudenko and Myshyakov (2009) 
was used in the present study: we take into account a full set of the evolu¬ 
tionary equations and allow the held variations not only in the inner volume 
but also on the lateral and top boundaries of the computational domain. 


6 
























distance from barycentre, arcsec 


Figure 3: Upper panel: dominant frequencies at three heights, and the magnetic field as a 
function of distance from the sunspot barycentre. Middle panel: dominant frequencies and 
the cut-off frequency derived from Equation 2. Bottom panel: magnetic field inclination 
angle derived from Equation 3 for the Hen 304A level, HMI Fel 6173A data, potential 
(B p ) and nonlinear force-free ( B np ) magnetic field extrapolation at the 2200 km level. 
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3. Results 


Figure 2 shows the spatial distribution of dominant frequencies in the 
sunspot at different heights from the photosphere to the corona. The black 
solid lines mark the inner and outer penumbra boundaries. Five-minute os¬ 
cillations concentrate in an annular area, that expands beyond the umbra 
boundaries with increasing height. Low-frequency oscillations appear in the 
lower photosphere of the penumbra. Such a topology resembles the con¬ 
cept explaining running penumbral waves as upwardly propagating waves 
that form a visual pattern (Rouppe van der Voort et ah, 2003; Bloomfield 
et ah, 2007; Kobanov et al., 2009). Low-frequency oscillations occupy the 
larger area in the picture plane at higher coronal level visible in Fe IX 171 A: 
the areas of dominant frequencies become elongated in the radial direction. 
Using the circular-shaped power distribution, one can plot frequency vs. dis¬ 
tance from the barycentre for different height levels observed with AIA and 
HMI. Figure 3 presents a dominant frequency [F(r)] averaged over the pixels 
located at distance r from the sunspot barycentre with a 0.6" step. The 
same averaging was used to obtain the longitudinal magnetic field [B,] and 
the absolute value of the magnetic-field inclination (Figure 3). We estimate 
the cut-off frequency f c (Figure 3, middle panel) and the magnetic-field in¬ 
clination angle [ip] (Figure 3, bottom panel) by Equations 2 and 3. Using 
potential B p and nonlinear force-free B np magnetic field extrapolation data 
we produce [<p\ vs. distance from barycentre plots for the 2200 km level. Both 
cases show a very similar dependence except for the outer penumbra (Fig¬ 
ure 3, bottom panel). In the Figure 3 middle panel one can see a penumbra 
region, in which the dominant frequency curve (304 A) follows the cut-off fre¬ 
quency curve. Earlier, Kuridze et al. (2009) showed that pressure impulses 
engendered in the photosphere can excite oscillations at the cut-off frequency. 
Consequently, the oscillations that coincide with the resonant frequencies of 
the cavities under the magnetic canopy amplify and can penetrate into the 
upper layers of the solar atmosphere (Srivastava et ah, 2008). 

High-frequency waves of the 6-7 mHz band concentrate at the sunspot 
central umbral part at all heights and form fragmented domains on the power 
maps (Figure 4). The power maps were produced by using entire time series. 
The NoRH at 17 GHz observation clearly show three-minute band oscillations 
in Stokes V polarisation channel (Figure 5). For the optically thin radiation 
the linc-of-sight magnetic field can be estimated according to Gelfreikh (2004) 
and is about 730 G. Extrapolation of the photospheric magnetic held to the 
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Figure 4: Intensity FFT power maps averaged in the 6-7 mHz band. The black solid lines 
mark the inner and outer penumbral boundaries. 
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Figure 5: Left panel: FFT power spectrum of the NoRH 17 GHz Stokes V signal. Right 
panel: FFT power spectrum of the Hen 304 A intensity. The signals were averaged over 
the area marked by the rectangle. The inner and outer penumbra boundaries are outlined 
by white solid lines. 



Figure 6: Intensity variations revealing wave trains in the three-minute band in the chro¬ 
mospheric and transition region lines and the NoRH 17 GHz Stokes V. 
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Figure 7: Fan-like structures in spatial localisations of low-frequency coronal oscillations 
(panels 1 and 2) and inverse sunspot image in the Fe IX 171 A band (panel 3). Panel 4: FFT 
power spectrum calculated over the whole field-of-view. The inner and outer penumbra 
boundaries are outlined by black and white solid lines. 
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Figure 8: Oscillation behaviour in coronal loops above sunspot. Left panel shows the 
intensity oscillation power spectra in the loop points marked with crosses in Figure 7. The 
results of frequency wavelet filtration are presented in the right panel. The thick lines 
show spectra and signals at the points located farther from the sunspot centre. 


upper atmosphere gives us the effective height of 2100-2200 km for the radio 
emission that corresponds to the transition region and the field strength of 
that value. 

Intensity variations in the narrow band of 3 minute for the 1600 and 
304 A bands show a wave train structure that is typical of sunspots. Similar 
structures are present in the NoRH 17 GHz polarisation signal (Figure 6). 
The measured time delay of three-minute propagating waves (Figure 6) from 
the level of 1000-1300 km (1600 A) to a height of 17 GHz emission gives a 
value of about 65 seconds, and 5 seconds from the level of 2200 km (304 A) to 
a height of 17 GHz emission. We can conclude that the calculated height of 
radio emission generation agrees with height estimations for other spectral 
lines (Avrett and Loeser, 2008; Zhang et ah, 1998) and with the analysis 
of characteristics of propagating three-minute waves (Reznikova et ah, 2012; 
Kobanov et ah, 2013). 

At the coronal heights, power maps reveal elongated structures that out¬ 
line the coronal loops best at the 1-2 mHz filtration band (Figure 7). Sim¬ 
ilarly to the two sunspots analysed in Kobanov et ah (2015), the 1-2 mHz 
band oscillations can be traced along the loops to find the time delay and 
propagating speed. For this purpose we find two points along the loop hav¬ 
ing almost identical FFT spectra. These points are marked with crosses in 
Figure 7. Then using the amplitude vs. time plot (Figure 6), we estimate 
the apparent propagation speed that appears to be ~45km/s. The angle be¬ 
tween the line of sight and the considered coronal fan was estimated from the 
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magnetic field extrapolation and was found to be about 40°. Taking into ac¬ 
count this value we find the deprojected wave propagation speed to be about 
70km/s. This value is lower than the expected sound speed in the corona. 
However, it is consistent with other measurements of the slow MHD wave 
propagation speed in coronal loops, ranging from 20km/s to 190km/s (See 
Table2 in De Moortel (2009)). It is interesting to note that low-frequency 
oscillations of 1-2 mHz band in sunspots can be connected with Evershed 
flow’s pulsations (Rimmele, 1995; Kobanov and Makarchik, 2004). The study 
of such a possible connection would require sunspots located near the solar 
disc limb. Other possible scenarios can be related to the stochastic motion 
of coronal loops or the non-damping transverse oscillations (Nistico et al., 
2013; Anfinogentov et ah, 2013). 

4. Conclusions 

The 5-7 mHz band high frequency oscillations concentrate inside small- 
scale elements in umbra and at the umbra-penumbra boundary. These kinds 
of oscillations show a well-known wave-train behaviour that also reveals it¬ 
self in the NoRH Stokes V 17 GHz data. The wave trains observed in UV 
and radio coincide with a cross-correlation coefficient up to 0.75-0.92. The 
time delays for the 1600 A, 304 A and the NoRH signal agree with the previ¬ 
ous results, obtained without radio data analysis, and indicate that the radio 
emission effective height is within the 2100-2200 km range which corresponds 
to the transition region and the magnetic field B~730G. This value agrees 
with one measured in He 110830 A line B~600-1500 G, (Solanki et ah, 2006) 
and is much higher than the estimated magnetic field at the coronal levels 
B~10G (Nakariakov and Ofman, 2001; Zaqarashvili et ah, 2013). The esti¬ 
mation of the magnetic fields at higher coronal levels would be possible with 
availability of new high cadence and frequency radio data from SSRT and 
CSRH instruments. 

The 1-2 mHz band low-frequency waves in NOAA 11711 reproduce the 
behaviour of NOAA 11479 and NOAA 11311 analysed by Kobanov et ah 
(2015): the power of these waves is concentrated in the annular zone around 
penumbra that expands with height. The estimated deprojected propagation 
speed of the wave train travelling along the single fan-loop is about 70km/s. 

The magnetic field inclination angle [p] in the sunspot was estimated in 
two independent ways: from the observations of the dominant oscillations 
frequency and by using the magnetic field extrapolation, [p] derived from 
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HMI photospheric data is ~35-40° for the umbra-penumbra boundary (12" 
from the sunspot centre), and corresponds to the held extrapolation results 
at the 2200 km level (Figure 3). [ 93 ] is obtained from indirect measurements 
for He II 304 A line is ~40-45°. The inclination estimation results fit in with 
the general structure of a sunspot and provide grounds for new research. 
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